#include "cHOdipole.h"

HOdipole::HOdipole(){
	setName("HO with dipole interaction");
	set_k(1.0);
	setUniverse(35);
	setWidth(getUniverse());
	setPotLeft(0.01);
	setPotRight(getUniverse());
        setUniLeft(0.01);
        setUniRight(getUniverse());
	setFx(getPotRight());
	setTx(0.01);
	setMinEnergy(1.5);
	//filter[activefilter]->calculate(this);
}

double HOdipole::getPotentialValueAt(double x){
	if(x>0) {
		return 0.5*get_k()*x*x + 2.0/(x*x);
	}
	else return 0;
}

HOdipole::~HOdipole(){
}
